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Abstract 

At the Large Hadron Collider (LHC) it will become possible for the first time to investigate 
experimentally the forward region in hadron-hadron collisions via high-p^ processes. In the LHC 
forward kinematics QCD logarithmic corrections in the hard transverse momentum and in the 
large rapidity interval may both be quantitatively significant. We analyze the hadroproduction 
of forward jets in the framework of QCD high-energy factorization, which allows one to resum 
consistently both kinds of corrections to higher orders in QCD perturbation theory. We compute 
the short-distance matrix elements needed to evaluate the factorization formula at fully exclusive 
level. We discuss numerically dynamical features of multi-gluon emission at large angle encoded in 
the factorizing high-energy amplitudes. 



1 



I. INTRODUCTION 



Experiments at the Large Hadron Collider (LHC) will explore the region of large rapidities 
both with general-purpose detectors and with dedicated instrumentation, including forward 
calorimeters and proton taggers [l|, 0, H, 0, H, @, 0]. The forward-physics program involves 
a wide range of topics, from new particle discovery processes [8, 0] to new aspects of 



strong interaction physics 0, OH l3] to heavy- ion collisions fl2l [l3 1 . Owing to the large 
center-of-mass energy and the unprecedented experimental coverage at large rapidities, it 
becomes possible for the first time to investigate the forward region with high-p^ probes. 

The hadroproduction of a forward jet associated with a hard final state X is pic- 
tured in Fig. [TJ The kinematics of the process is characterized by the large ratio of sub- 
energies si/s ^> 1 and highly asymmetric longitudinal momenta in the partonic initial state, 
k\ ■ P2 ^ &2 " Pi- At the LHC the use of forward calorimeters allows one to measure events 
where jet transverse momenta px > 20 GeV are produced several units of rapidity apart, 
Ay > 4 -T- 6 [H, H, 0] . Working at polar angles that are small but sufficiently far from the 
beam axis not to be affected by beam remnants, one measures azimuthal plane correlations 
between high-p^ events widely separated in rapidity (Fig. [2]). 
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FIG. 1: Jet production in the forward rapidity region in hadron-hadron collisions. 



The presence of multiple large- momentum scales implies that, as was recognized 



in 



141 . Il5l . [16| ], reliable theoretical predictions for forward jets can only be obtained after 



summing logarithmic QCD corrections at high energy to all orders in a s . This has motivated 
efforts [HHH, ]20J to construct new algorithms for Monte Carlo event generators capable 
of describing jet production beyond the central rapidity region. Note that an analogous 
observation applies to forward jets associated to deeply inelastic scattering [2l|, |22|. Indeed, 
measurements of forward jet cross sections at HERA j23[ indicate that neither fixed-order 
next-to-leading calculations nor standard shower Monte Carlo generators [TJ 



23. 24 



e.g. 

Pythia or Herwig, are able to describe forward jet ep data. Improved methods to eval- 
uate QCD predictions are needed to treat the multi-scale region implied by the forward 
kinematics. 

In this work we move on from the observation that realistic phenomenology in the LHC 
forward region will require taking into account at higher order both logarithmic corrections 
in the large rapidity interval (of high-energy type) and logarithmic corrections in the hard 
transverse momentum (of collinear type). The theoretical framework to resum consistently 
both kinds of logarithmic corrections in QCD perturbation theory is based on high-energy 
factorization at fixed transverse momentum [251 ] . 

This formulation depends on unintegrated distributions for parton splitting, obeying ap- 
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FIG. 2: (Left) High-p^ events in the forward and central detectors; (right) azimuthal plane seg- 
mentation. 



propriate evolution equations, and short-distance, process-dependent matrix elements. The 
unintegrated-level evolution is given by evolution equations in rapidity, or angle, param- 
eters. Different forms of the evolution, valid in different kinematic regions, are available. 



See 0, H HQ, and references therein, for recent work in this area and reviews. The 
short- distance matrix elements, needed in the evaluation of the factorization formula, are 
the subject of this paper. We obtain their explicit expressions in a fully exclusive form, in- 
cluding all partonic channels, and present results of numerically integrating them over final 
states. Such matrix elements, though not on shell, are gauge invariant and perturbatively 
calculable. They factorize in the high energy limit in front of (unintegrated) distributions 
for parton splitting not only in the collinear emission region but also at finite angle. In par- 
ticular, they can serve to take into account effects of coherence from multi-gluon emission, 
away from small angles, which become important for correlations among jets across long 
separations in rapidity. We give a numerical illustration of the high-A;^ behavior resulting 
from such finite-angle radiation. 

On one hand, once convoluted with the small-x gluon Green's function according to the 



method |25l . |30|, these matrix elements control the summation of high-energy logarithmic 



corrections, contributing both to the next-to-leading-order BFKL kernel 3JJ and to the jet 



impact factors [32|, [33J]. On the other hand, they can be used in a shower Monte Carlo 



implementing parton-branching kernels at unintegrated level (see e.g. |34l . |35| for recent 
works) to generate fully exclusive events. We leave these applications to a separate paper. 

The paper is organized as follows. After recalling the factorized form of the cross sections 
in Sec. HIl we present the high-energy amplitudes in Sec. 1111} and discuss basic properties 
and numerical results in Sec. [IVl We summarize in Sec. [V] 



II. HIGH-ENERGY FACTORIZED CROSS SECTIONS 



High-energy factorization [25| allows one to decompose the cross section for the process 
of Fig. [Tj into partonic distributions (in general, unintegrated) and hard-scattering kernels, 
obtained via the high-energy projectors |25l. l30f] from the amplitudes for the process P1+P2 —> 
P3 + Pi + 2 massless partons. The basic structure is depicted in Fig. [3j 

With reference to the notation of Fig. [3], let us work in the center of mass frame of the 
incoming momenta 



Pi 



V^/2(l,0,0r) 



p 2 = y/3/2{0 i l,0 T ) 



2pi • P2 = S 



where, for any four- vector, p^ = (p + ,p~ ,pr), with p^ = (p°±p 3 ) / y2 and px two-dimensional 
euclidean vector. Let us parameterize the exchanged momenta in terms of purely transverse 
four-vectors k± and k±i and longitudinal momentum fractions & and ^ as 



Pi~P5 = h= £m + + £i£>2 , P2 - Pe = h = &P2 + k± + i 2 px . 



(2) 
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For high energies we can introduce strong ordering in the longitudinal momenta, £1 ^> |£ 2 |> 
£2 3> Further, we make the forward region approximations (p4 + P6) 2 ^ (P3 +P4) 2 , 
kx ~ ^ipi, A; 2 ~ + so that 

P5 ^ (1 - 6)pi , Pe ^ (1 - 6)P2 - fcj. , 6>6- (3) 

It is convenient to define the rapidity-weighted average of dijet transverse momenta, 

Q T = (1 - u )p T4 - upT3 ; where v = (p 2 Pi) / [(P2 Pi) ~ (P2P5)] , ( 4 ) 
and the azimuthal angle 

cosy? = Q T ■ k T /\Q T \\k T \ ■ (5) 
We consider the differential jet cross section in Qt and (p. 




(a) 



FIG. 3: (a) Factorized structure of the cross section; (b) a graph contributing to the qg channel 
matrix element. 



According to the factorization [25|, [32], the jet cross section can be computed as (Fig. [3^) 



= £ J dix d& d 2 k T cf> a/A (ti) fcr, Qr, fcr) , (6) 

where the sum goes over parton species, <fi are the parton distributions defined from the 



unintegrated Green's functions introduced in [30| for both gluon and quark cases, and a 
is the hard cross section, calculable from the high-energy limit of perturbative amplitudes 
(Fig.Eb). 

The physical picture underlying Eq. ([6]) is based on the fact that initial-state parton 
configurations contributing to forward production are asymmetric, with the parton in the 
top subgraph being probed near the mass shell and large x, while the parton in the bottom 
subgraph is off-shell and small-x. Eq. ([6]) embodies this picture through the longitudinal 
and transverse momentum dependences of both and a. 

For phenomenological studies we will be interested in coupling Eq. ([6]) to parton showers 
to achieve a full description of the associated final states. To this end we need the matrix 
elements defining the hard-scattering kernels in a fully exclusive form. We give the results 
in the next section. 
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III. MATRIX ELEMENTS FOR FULLY EXCLUSIVE EVENTS 



The matrix elements determining the hard-scattering kernels a can be viewed as a suitably 



defined off-shell continuation of scattering amplitudes at lower order |25j. They can be 
obtained by applying to scattering amplitudes M. the high-energy eikonal projectors [25|, 30], 

M H = pH^... Mm (ku ^ {pi}) ) pH^... K (?) 

V fc ±l fc _L2 

Although they are not evaluated on shell, they are gauge invariant and their expressions are 
simple. The utility of these matrix elements is that in the high-energy limit they factorize 
not only in the collinear emission region but also in the large-angle emission region. As 
long as the factorization is carried out in terms of distributions for parton splitting at fixed 



transverse momentum, as in Eq. they can be useful to include coherence effects [32 
from multi-gluon emission across large rapidity intervals, not associated with small angles. 
The results for the matrix elements in exclusive form are given by 

M qg ^ qg = C 1 Af> + C 1 Wr b) , M gg ^ q - q = C 2 A 2 ab) +C 2 A [ 2 nab) , M„ = C 3 A 3 (8) 



where 



A (ab) = I hk 2 \ (k lP2 ) 2 + {p 2 p 3 f , g x 
1 hp 2 (&1P4) O3P4) 



A (nab) = f h h \ (hp 2 ) 2 + (>2P3) 2 / (P3P4) (hP2) + (hpi) (P2P3) _ A ^ 
\klP2) (fcip 4 ) {PzPi) \{hP3)(P2P4) {hP3){P2Pi) J 

,(ab) _ (h h\ (P2P3) 2 + (P2P4Y nn 
2 " \hp 2 ) (k lP4 ) {k lP ,) ' 

A (nab) = ( h k 2 \ (P2P3) 2 + (P2P4) 2 ( (fcip 4 ) {P2P3) jklPz) (P2Pi) _ A < l2 , 
\klP2J (hpi) (hp 3 ) V G°3P4) {kip 2 ) {P3Pi){k 1 p 2 ) ) 

. fk 1 k 2 \ (P3P4)(k 1 p 2 ) + (k 1 p i )(p 2 p 3 ) + (p 2 p 4 )(kip 3 ) r 4 4 4 , 

^ = (jW (, 2 ,4)(^P4)(P3,4)(^, 2 )(, 2 P3)(^,3) ^ + ^ ] ' 

(13) 

and Q = g\N 2 c - 1)/(4A C 2 ), C 1 = C l C A /{2C F ), C 2 = g A /(2N c ), C 2 = C 2 C A /(2C F ), 
C 3 = g*N 2 /(N*-l). 

The results above contain the dependence on the transverse momentum kj_ along the 
parton lines that connect the hard scatter to the parton distributions. Nevertheless, they 
are short- distance in the sense that they can be safely integrated down to k_L = 0. That 
is, the high-energy projection is designed so that all infrared contributions are factored out 
in the nonperturbative Green's functions <fi in Eq. ([6]). An explicit numerical illustration is 
given in the next section. 1 



1 Although the hard-scattering functions constructed from the amplitudes in Eqs. ((9|l- (fT3|) are not coefficient 
functions in the conventional sense of the operator product expansion, they can be related to such objects, 



functions in the sense of the high-energy OPE of 36 1. 



for inclusive variables, along the lines e.g. of 30 L They could be interpreted in terms of coefficient 
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The role of Eqs. (l9l)- (jl~3l) is twofold. On one hand, they give the high-energy limit of multi- 
parton matrix elements in the forward region, which may be of direct phenomenological 



significance. On the other hand, because of the factorization theorem |25|], logarithmically 
enhanced corrections for large rapidity can be systematically obtained to all orders in a s 
from those in the (unintegrated) distributions for parton splitting once the hard scattering 
functions are known at finite kj_. To this end the detailed form of the fall-off at large |kj_ | 2 
is relevant. 

In the next section we discuss the behavior at high transverse momentum numerically. 
This behavior reflects properties of gluon emission at large angle encoded in the high-energy 
amplitudes. These are relevant, along with large- angle effects in the Sudakov region (see 
e.g. [3~7|). to achieve a full treatment of gluon coherence effects 38] capable of describing jet 
final states across the whole rapidity phase space. A uniform treatment of the high-energy 
and Sudakov regions is still an open issue [39| of interest for parton-shower implementations. 



IV. NUMERICAL RESULTS 

We now partially integrate the amplitudes over final states. We work at the level of 
hard scattering matrix elements, leaving the treatment of parton evolution by showering to 
a separate study jioj. We concentrate on the region of hard emissions, where jets are well 
separated. Regions near the boundary of the angular phase space are sensitive to infrared 
radiation and can be addressed within a full parton-shower description of the process. 

We consider the differential distribution in the transverse variable Qt and azimuthal 
angle tp. The variable Qt describes the imbalance in transverse momentum between the 
hardest jets, weighted by u, according to Eq. In Figs. H] and Owe show numerical results 
versus transverse momentum and versus energy (qg channel). 




FIG. 4: The fcr/Qr dependence of the factorizing qg hard cross section at high energy: (left) C\ 
term; (right) C F C A term (Ci^-S'/Qt = 1q2 > a s = 0-2). 

The curves in Fig. H] measure the kx distribution of the jet system recoiling against the 
leading di-jets. The result at kr/Qr ^ in these plots returns the lowest-order result, i.e., 
the leading-order process with two back-to-back jets, 

d^^ a "f (0) (PT/ s ) > Qt - Pt = \Pts\ = M , (14) 
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where s = (p 3 + P4) 2 , and /( ) is given by 



/<">(*) = _L_ [C&*(l + z) + 2C F C A (1 -3z + z 2 )] . (15) 

The dependence on kx and <p plotted in Fig. H] is the result of higher-order gluon radiation, 
treated according to the high-energy asymptotics. The different behaviors in cp for the C F 
and CaCf terms reflect the fact that the former comes from the insertion of gluons on 
fermion-exchange amplitude while the latter comes from the insertion of gluons on vector- 
exchange amplitude. 

Fig. [5]shows the energy dependence for fixed kp/Qr- The constant asymptotic behavior at 
large s due to color-octet spin-1 exchange distinguishes the CfCa term from the Cp term. 
The dependence on the azimuthal angle in Figs. 0] and [5] is relevant, especially because 
forward jet measurements will rely on azimuthal plane correlations between jets far apart in 
rapidity (Figj2j). 




FIG. 5: The energy dependence of the qg hard cross section {Ht/Qt — !)• 



While Eq. (fTST) gives the collinear emission limit, we see from Fig. H] that multi-gluon 
radiation at finite angles sets a dynamical cut-off at values of kr of order Qj>, 



k T < /i 



c Qq 



(16) 



The physical meaning of this result is that the summation of the higher-order logarithmic 
corrections for large y ~ In s/p T is precisely determined (25I . 32| by convoluting the uninte- 
grated splitting functions over the kp- dependence in Fig. HI via the distributional relation 



d 2 kn 



A 7 1 



J P(^t) - @(/i - k T ] 



(17) 



So the results in Fig. H] illustrate quantitatively the significance of contributions with k? — 
Qt in the large-?/ region. Non-negligible effects arise at high energy from the finite-Zcr 
tail. These effects are not included in collinear-branching generators (and only partially in 
fixed-order perturbative calculations), and become more and more important as the jets are 
observed at large rapidity separations. 

Observe that calculations based on the unintegrated formalism will in general depend 
on two scales, \x and the rapidity, or angle, cut-off 0, [T^ 



41] . See e.g. the one-loop 
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calculation 29 



for an explicit example. It will be of interest to investigate the effect of 
Eq. (fT6|) on the behavior in rapidity distributions 4^ . 

Results for gluon-gluon channels are reported in Fig. [61 Note the large effect of the purely 
gluonic component. The behavior versus hr is qualitatively similar to that in Fig. [H Calcu- 
lations in progress jio|, including parton showering, indicate that quark and gluon channels 



give contributions of comparable size in the LHC forward kinematics. The inclusion of both 



is relevant for realistic studies of phenomenology 11, 42| . Since the forward kinematics 
selects^ asymmetric parton momentum fractions, effects due to the x — > 1 endpoint behav- 
ior 
well 



29l . l43j at fixed transverse momentum may become phenomenologically significant as 




FIG. 6: The kx/Qr dependence of the factorizing gg matrix elements: (a) gg — ► gg; (b) gg — > qq 
Cf term; (c) gg — > qq Ca term (Zi&S/Qt = 10 2 , a s = 0.2). 

We conclude this section by recalling that dynamical effects of high parton densities have 
been studied 0, |44j as potential contributions to forward jet events. We note that if such 



effects show up at the LHC, the unintegrated formulation discussed above would likely be 
the natural framework to implement this dynamics at parton-shower level. 



V. CONCLUSIONS 

Forward + central detectors at the LHC allow jet correlations to be measured across 
rapidity intervals of several units, Ay > 4 -j- 6. Such multi-jet states can be relevant to new 
particle discovery processes as well as new aspects of standard model physics. 

Existing sets of forward-jet data in ep collisions, much more limited than the potential 
LHC yield, indicate that neither conventional parton-showering Monte Carlo generators nor 
next-to- leading-order QCD calculations are capable of describing forward jet phenomenology. 
These observations motivate studies of improved methods to compute QCD predictions in 
the multiple-scale kinematics implied by the forward region. 

We have analyzed the high-energy factorization that serves to sum consistently to higher 
orders in a s both the logarithmic corrections in the large rapidity interval and those in the 
hard jet transverse energy. We have determined the gauge-invariant (though not on shell) 
high-energy amplitudes, which are needed to evaluate the factorization formula for forward 
jet hadroproduction. 

Our results can be used along with k^-dependent kernels for parton branching. They can 
serve to construct predictions for exclusive observables associated to forward jets, including 
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jet correlations, that take into account gluon coherence not only in the collinear emission 
region but also in the large-angle emission region. 
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